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The spin axes of gyroscopes experimentally define local non-rotating frames, i.e. the time- 
evolution of axes of inertial frames. But what physical cause governs the time-evolution of gyroscope 
axes? We consider linear perturbations of Friedmann-Robertson- Walker (FRW) cosmologies with 
k = 0, i.e. spatially flat. We ask: Will cosmological vector perturbations (i.e. vorticity or rotational 
perturbations) exactly drag the spin axes of gyroscopes relative to the directions of geodesies to 
quasars in the asymptotic unperturbed FRW space? Using Cartan's formalism with local orthonor- 
mal bases we cast the laws of linear cosmological gravitomagnetism into a form showing the close 
correspondence with the laws of ordinary magnetism. Our results, valid for any equation of state 
and any form of the energy-momentum tensor for cosmological matter, are: 1) The dragging of a 
gyroscope axis by rotational perturbations of matter beyond the //-dot radius from the gyroscope 
is exponentially suppressed, where H is the Hubble rate, and dot is the derivative with respect to 
cosmic time. 2) If the perturbation of matter is a homogeneous rotation inside some radius around a 
gyroscope, then exact dragging of the gyroscope axis by the rotational perturbation is reached expo- 
nentially fast as the rotation radius gets larger than the H-dot radius. 3) For the most general linear 
cosmological perturbations the time-evolution of all gyroscope spin axes and the axis directions of 
all local inertial frames exactly follow a weighted average of the rotational motion of cosmological 
matter, i.e. there is exact frame-dragging everywhere. The weight function is the density of mea- 
sured angular momentum of matter times (1/r) times the Yukawa force (— d/dr)[(l/r) exp(— fir)], 
where r is the geodesic distance from the source to the gyroscope. The exponential cutoff is given by 
(i 2 — —4(dH/dt). Except for the Yukawa cutoff the weight function is the same as in the integrated 
form of Ampere's law. — Our results demonstrate (in first-order perturbations of any type for FRW 
cosmologies with k = 0) the validity of Mach's hypothesis that axes of local non-rotating frames 
precisely follow an average of the motion of cosmic matter. 

PACS numbers: 04.20.Cv, 04.25.Nx, 98.80.Jk 



I. INTRODUCTION 

A. The observational fact 

In tests of general relativity in the solar system two 
types of things are compared: 1) On the one hand mea- 
surements of the precession of perihelia or of gyroscopes' 
spin axes in Gravity Probe B [l| relative to distant stars 
and quasars. 2) On the other hand the predictions from 
the solutions of Einstein's equations for the solar system 
in asymptotic Minkowski space, which do not contain dis- 
tant stars or quasars. In this comparison the assumption 
is made (often not spelled out) that in the solution of 
Einstein's equations for the solar system the asymptotic 
Minkowski space is nonrotating relative to distant stars 
and quasars. This implicit assumption is tested to high 
accuracy by the comparison of the observed perihelion 
precession with the prediction of general relativity. This 
implicit assumption is a basic observational fact, which 
has been called "Mach 0" e.g. in J|. — Distant stars and 
quasars have proper motions relative to the uniform Hub- 
ble flow, but the angular parts are unmeasurable today. — 
We conclude that the measurements of perihelion shifts 
and the measurement undertaken by Gravity Probe B 
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are tests of two things combined, on the one hand tests 
of Einstein's equations, on the other hand tests of the 
principle "Mach 0" . 

No explanation for the observational fact "Mach 0" 
is given in classical mechanics, special relativity, and 
general relativity for isolated systems in asymptotic 
Minkowski space, except by invoking an accident of initial 
conditions. In these theories one could have different ini- 
tial conditions, where distant stars and quasars could be 
orbiting around us relative to our gyroscopes. — Within 
these three theories (and in the generic case of general 
relativity) the local non-rotating frames (axis directions 
of inertial frames) can be experimentally determined by 
axes of gyroscopes. Conversely the time evolution of gy- 
roscope axes is dictated by the laws of inertia, i.e. the 
gyroscope axes cannot rotate with respect to local iner- 
tial axes. This is self-consistent: One uses two gyroscopes 
(with spin axes not aligned) to experimentally construct, 
i.e. operationally define, the local non-rotating frame, 
and one uses all other local gyroscopes and other exper- 
iments (e.g. Foucault's pendulum) to test the laws of 
inertia with respect to rotations. — But the question re- 
mains: What physical cause governs the time-evolution 
of the axes of inertial frames and of gyroscopes? These 
three theories do not contain such a physical cause, quite 
to the contrary they contain an absolute element, namely 
the non-rotating frame in classical mechanics and spe- 
cial relativity, and the asymptotic non-rotating frame 
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in general relativity for isolated systems in asymptotic 
Minkowski space. This absolute element is experimen- 
tally established by Newton's bucket experiment and to- 
day by gyroscopes or the Sagnac effect in inertial guid- 
ance systems. This absolute element is equivalent (with 
respect to rotations) to Newton's absolute space. 

B. Mach's Principle 

In the 1870's and 1880's Ernst Mach Q 0] stated force- 
fully the hypothesis (as an alternative to Newton's abso- 
lute space) that the axes of non-rotating frames (i.e. axes 
of gyroscopes) in their time-evolution are determined by, 
are exactly dragged by, precisely follow "some average" of 
the motion of matter in the universe. Mach's formula- 
tion is given in yj. See also @ — Today one generalizes 
"matter" to "matter-energy", i.e. the principle states 
that there is exact frame-dragging for the axis directions 
of inertial frames by "some average" of the energy cur- 
rents in the universe. This is what we take as the formu- 
lation of Mach's principle at the level of linear cosmolog- 
ical perturbations. At the level of linear perturbations 
gravitational waves do not carry energy in the sense of a 
pseudo-energy-momentum tensor. — Mach's Principle in 
this formulation follows from general relativity for linear 
perturbations of FRW universes with k — 0, as will be 
demonstrated by our result Eq. J3J. 

Many alternative versions have been proposed under 
the name of Mach's principle by other authors after Mach 
0,0- We shall di scuss Einstein's proposal of 1918 Q 
at the end of sect. IXHI An incomplete list of ten in- 
equivalent versions has been given and briefly discussed 
in 0- Six of these versions are definitely not satisfied 
for Einstein gravity in a cosmological context, and these 
six versions have absolutely no counterpart in the writ- 
ings of Mach. Another version is satisfied in Einstein 
gravity, even in a non-cosmological context, but it is too 
weak: partial dragging as opposed to exact dragging, i.e. 
inertial frames influenced as opposed to fully determined 
by the motion of cosmic matter. Because of the many 
conflicting versions proposed under the name of Mach's 
principle by other authors after Mach, some physicists 
have been under the impression that the issue is "con- 
fusing and ill-defined". — The version, "the theory con- 
tains no absolute elements," explained by J. Ehlers in 
[Toj , is intimately connected to Mach's starting point (no 
absolute space) and formulated in a modern way. This 
version is a general requirement, but it is not a specific 
implementation in contrast to our result Eq. J2J, which 
gives an explicit demonstration of how Mach's principle 
is implemented within general relativity. 



C. Gravitomagnetism 

At the time of Mach there was no known mechanism 
by which matter in the universe could influence the mo- 



tion of gyroscope axes. With General Relativity came the 
needed mechanism, gravitomagnetism. Thirring in 1918 
[Tl| considered a rotating infinitely thin spherical shell 
with uniform surface mass density and total mass M, and 
he analyzed the partial dragging of the axes of inertial 
frames inside the rotating shell for points near the ori- 
gin. In the weak field approximation, Gn(M/ -R) s hcii 1, 
and to first order in the angular velocity f2 s heii he found 
that the axes of local inertial systems (and gyroscopes) 
at all points near the origin rotate relative to asymp- 
totic Minkowski space with the same angular velocity 
^gyro = /drag^sheii- For the dragging fraction / drag he 
obtained /drag = ^G^(A1/R) shc ii <C 1. — Gravitomag- 
netic effects in the solar system are exceedingly small. 
The Lense- Thirring effect by the rotating earth on the 
spin axes of gyroscopes on Gravity Probe B is predicted 
to be only 43 milli-arc-sec per year — Brill and Co- 
hen ^3 and Lindblom and Brill ^| analyzed dragging 
by spherical shells to first order in the angular velocity 
but to all orders in Gn(M/R) shell- They found that in- 
side the shell one has a Minkowski space which rotates 
relative to asymptotic Minkowski space. In the limit of 
(M/R) s heii approaching the value for a black hole they 
found perfect dragging of inertial axes inside the rotating 
shell. 

Exact dragging by the masses in the universe is re- 
quired by Mach's principle, not a small influence. There- 
fore one must go to cosmological models. A simple, cos- 
mologically relevant model is a uniform rotational mo- 
tion of cold matter with f2 ma ttcr and mass density p 
both constant out to a radius i? ro t around a gyroscope 
and f2 m atter = outside R-cot- Starting with a weak 
field perturbation on a Minkowski background this simple 
model is a superposition of a sequence of Thirring shells. 
This gives a dragging fraction which grows quadratically 
in R rot , /drag = 2G N (M/R) rot = (8n /3)G N pR 2 ot . The 
Friedmann equation for k = reads H 2 — (8ir/3)Gp, 
hence 2G^Mh/Rh = 1, where Mjj is the mass inside 
the Hubble radius Rh- This gives /drag = R 2 ot /R 2 H - Al- 
though the weak field perturbation on a Minkowski back- 
ground implies /drag = 2G N (M/R) rot < 1, applying the 
formula beyond its region of validity would give /drag — * 1 
for i? rot — > Rh, a "satisfying degree of self-consistency" 
in the words of Misner, Thorne, and Wheeler 0], and 
/drag 2> 1 for R rot ^> R H , the problem of 'overdragging'. 
What is needed is an analysis using cosmological per- 
turbation theory including superhorizon scales. This is 
given in this paper. 



II. SUMMARY AND CONCLUSIONS 

In this paper we analyze realistic cosmological mod- 
els. This is in contrast to toy models involving one or 
many shells (embedded in FRW models), infinitely thin 
and rotating around axes through one point, which are 
contained in [r|, [r|, [T^J. In models [ljl, Lu| "almost 
all the matter of the universe" must be "redistributed on 
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the shell" in order to make frame dragging almost per- 
fect, i.e. these toy models are extremely far from phys- 
ical cosmology. Also we consider realistic cosmological 
matter of any type as opposed to the contrived energy- 
momentum tensors discussed in the literature. We start 
from a Friedmann-Robertson- Walker (FRW) cosmology 
which is spatially flat (i.e. k = 0). We first add the most 
general linear cosmological perturbations in the vorticity 
sector (vector perturbations), i.e. perturbations derived 
from divergenceless vector fields. We ask: Will rotational 
motions of cosmological matter exactly drag the axes of 
any gyroscope relative to the directions of geodesies from 
the gyroscope to quasars in the asymptotic unperturbed 
FRW space? In Sec. XII we also include scalar and tensor 
perturbations. 

A conference paper about our work appeared in [Tsf 
and a preliminary version of a full article in 19] . A re- 
action 20] to our conference paper will be discussed in 
sect- IVim 

Our first result is that the dragging of axes of gyro- 
scopes (i.e. axes of local inertial systems) by matter 
beyond the -ff-dot radius RdH/dt is exponentially sup- 
pressed, where RdH/dt = {—dH/dt)^ 1 / 2 , H is the Hubble 
rate, and dot is the derivative with respect to cosmic 
time. In this sense physics here is decoupled from the 
asymptotic FRW universe. 

Our second result refers to the simple model discussed 
above (uniform rotation of matter inside a radius i? ro t 
around a gyroscope) . It states that the dragging fraction 
/drag approaches the value 1, i.e. exact dragging, expo- 
nentially fast as i? ro t gets larger than the H-dot radius. 
The dragging fraction as a function of R rot is shown in 
Fig. 1. The problem of 'overdragging' discussed in 0] is 

f [ drag ] 

/ 




R [ rot ] 
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FIG. 1: The dragging fraction /d rag ^gyro — /"drag ^matter 
versus i? ro t for a uniform rotation of cosmological matter in- 
side a radius R TO t around a gyroscope. The solid curve shows 
the result of cosmological perturbation theory: As R rot grows 
beyond the .ff-dot radius, /drag approaches the value 1 ex- 
ponentially fast, i.e. exact dragging. The dashed curve is 
based on extrapolating the results of gravitomagnetic pertur- 
bation theory on a Minkowski background beyond the region 
of validity: It shows the problem of 'overdragging', which is 
removed by the exponential suppression of contributions from 
super- H-dot scales. 



removed by the exponential suppression at super- H-dot 
scales, which is missing in |16) and other papers. — For 
linear cosmological perturbations both £! ma tter and r2 gyr o 
are infinitesimally small, but their ratio approaches the 
value 1 for i? rot much larger than the H-dot radius. 

Our third result concerns the most general vorticity 
perturbation in linear approximation. To understand the 
result one needs two crucial inputs from cosmological per- 
turbation theory in the vector sector (see Sec. Ill): 

1. The slicing of space-time in slices St in the vorticity 
sector is unique, as emphasized by Bardeen plj . 

2. The intrinsic geometry of 3-space remains unper- 
turbed, i.e. for vorticity perturbations of FRW with 
k = each slice S t remains a Euclidean 3-space 
with its global parallelism. 

Our result states what specific average of the energy flow 
J e in the universe determines the precession rate O gyr o of 
gyroscope axes, 

^gyro - ~B s {r = 0) = 

2G N /dV{I[r'xJ e (f')]}{^^-}, (1) 

where r ' is the position of the source relative to the gy- 
roscope. J e is the energy current, with J e = {p + p)il 
for perfect fluids, and B g is the gravitomagnetic field, of 
which the general operational definition is given in sec- 
tion llVI All quantities in Eq. Q are directly measurable, 
i.e. they have a gauge invariant meaning, as is explained 
after Eq. (|42|l . — The weight function is the density of 
measured angular momentum of matter times (1/r) times 
the Yukawa force T M (r) = (— d/dr)[(\/r) exp(— fir)]. 
The exponential suppression is determined by fi 2 = 
—A{dH/dt) — 167tGn(p + p)- Except for the exponential 
factor the weight function is the same as in the integrated 
form of Ampere's law for ordinary magnetism. 

Note the fundamental difference between cosmological 
gravitomagnetism, Eq. and Ampere's law. The latter 
does not hold in a rotating reference frame, unless one 
introduces fictitious forces. In contrast Eq. remains 
valid, as it stands, in a frame which is rotating with an 
angular velocity Q* relative to asymptotic quasars. This 
holds because both sides of Eq. Q change by the same 
amount: On the right-hand side, (r' x J c ) changes by 
(p + p)[f ' x (£1* x r')]. Carrying out the integration gives 
(l* with a prefactor 1 for the change on the right-hand 
side, which is equal to the change on the left-hand side. — 
The fact that Eq. JQl holds in any frame which is rotat- 
ing relative to asymptotic quasars establishes that the 
asymptotic inertial frame has no influence in cosmologi- 
cal gravitomagnetism; the time-evolution of local inertial 
axes is determined exclusively by the weighted average of 
cosmological matter flows. 

Eq. Q can be rewritten to show explicitely that the 
precession of any gyroscope, ^ gy ro, exactly follows, i.e. 
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is exactly dragged by a weighted average of fi ma tter with 
a weight function W{rpQ), where P and Q are the po- 
sitions of the gyroscope and the source. For reasons of 
symmetry under rotations and space reflection, the ve- 
locity field on a shell of a given radius rpQ can only con- 
tribute to the gyroscope's precession through its term 
with {£ = 1, odd parity sequence}, which is equivalent 

to a rigid rotation with the angular velocity f2 m atter(?")- 
We obtain 

/>oo 

S\yro = dr fimattcrW W{r), 

JO 

W(r) = i/i 2 r 3 F M (r) (2) 

This is our most important result. It shows that ^ gyr o is 
the weighted average of f2 m attcr, i-e. the evolution of in- 
crtial axes exactly follows the weighted average of cosmic 
matter motion. The weight function W(r) has its nor- 
malization (integral over r) equal to unity, as it must be 
for an averaging weight function in any problem. That 
the integral of the weight function is equal to unity cru- 
cially depends on the exponential cutoff in the Yukawa 
force Y^r). Details about Eqs. J2J are given in the para- 
graph with Eqs. ijTfl) and (|4*5)l . — Our results (valid for 
linear perturbations of a FRW universe with k = and 
for any equation of state), particularly Eq. J2J), are a clear 
demonstration of how Mach's principle is implemented 
within cosmological general relativity. 

Our aim is to obtain the laws of linearized cosmolog- 
ical gravitomagnetism in a form which shows the corre- 
spondence with electromagnetism in a (3+l)-dimensional 
split. In sect. IV we give the general operational defini- 
tion for the gravitoelectric field E g and the gravitomag- 
netic field B g via measurements by any choice of fidu- 
cial observers (FIDOs). These operational definitions are 
generally valid, i.e. beyond perturbation theory. The op- 
erational definitions need Cartan's formalism with local 
orthonormal bases (LONBs) and the corresponding fidu- 
cial observers. For our specific problem we choose FIDOs 
with spatial basis vectors fixed to directions of geodesies 
on S t from the FIDO to quasars in the asymptotically 
unperturbed FRW universe. This construction is clearly 
independent of the chosen coordinate system, i.e. it is a 
gauge- invariant construction. 

In sect. V we first give the general method to com- 
pute connection coefficients for local orthonormal bases 
in Cartan's formalism. We apply this general method 
to obtain the connection coefficients for linear cosmolog- 
ical gravitomagnetism and the equations of motion for 
free-falling test particles. 

In sect. VI we compute curvature in Cartan's for- 
malism. The crucial equation for cosmological gravito- 
magnetism is Einstein's Ggj-equation for vorticity per- 
turbations, the momentum constraint. It has the form 
of Ampere's law, except that for vorticity perturbations 
of FRW space there is an additional term — 4(dH/dt)A g , 
which is responsible for the exponential suppression of 



super-77-dot scales, 

curing - A{dH/dt)A g = -16ttG n J s . (3) 

In contrast to the Ampere-Maxwell equation of ordinary 
electrodynamics, the Maxwell term d t E is absent in the 
time- dependent context of gravitomagnetodynamics; see 
the comments after Eq. (|36[) . The solution of Eq. © is 
given in Eq. (Q. 

In sect. VII we give a detailed discussion of our 
three main results for Mach's principle, which were sum- 
marized above. — We then discuss Einstein's objection: 
"Mach conjectured that inertia would have to depend 
upon the interaction of masses ... and their interactions 
as the original concepts. The attempt at such a solution 
does not fit into a consistent field theory." We point out 
the remarkable fact that the solution of Einstein's local 
field equation for Ggj, i.e. of the momentum constraint 
equation J2Jl, has the form of an instantaneous action-at- 
a-distance, Eq. Q). See also refs. p|, |2^ . 

In sect. VIII we explain that the density of measured 
angular momentum, formed from the LONB components 
T?, is the relevant input for Mach's principle on the right- 
hand side of Einstein's equation, and not the density of 
conserved angular momentum, e.g. the coordinate-basis 
component T^, which was considered by other authors, 
e.g. ^(|. This is the reason why these and other authors 
did not obtain the Yukawa suppression beyond the iJ-dot 
radius. This Yukawa suppression is crucial for obtaining 
the normalization to unity in the weight function W(r) 
in Eq. (j2J), and it is crucial to remove the problem of 
'overdragging'. 

In sect. IX we discuss Einstein's objection "If you have 
a tensor T^ v and not a metric, then this does not mean- 
ingfully describe matter The statement that matter 

by itself determines the metric ... is meaningless." From 
this objection Einstein drew the conclusion that "one 
should no longer speak of Mach's principle." — His objec- 
tion applies to coordinate-basis components T^ v , but we 
explain why it does not apply to LONB components T.t. 
LONB components are determined directly by matter 
measurements using only rj^ — diag{ — 1, +1, +1, +1 }, 
which is available before having solved Einstein's equa- 
tions. In contrast T M „ needs matter measurements plus 
g^u, which is available only after having solved Einstein's 
equations. 

In sect. X we explain that the vanishing of local vor- 
ticity measured by non-rotating observers is not relevant 
as a test of Mach's principle, because Mach's principle 
requires an average over all matter in the universe with 
a given weight function. Therefore we conclude that the 
claim by Ozsvath and Schiicking that their Bianchi IX 
model violates Mach's principle is not conclusive. 

In sect. XI we derive the boundary conditions needed 
when solving Einstein's equations for a system in asymp- 
totic Minkowski geometry, if one does not restrict the 
allowed coordinate systems to those which are asymptot- 
ically non-rotating. The boundary conditions are needed 
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to encode the effects of sources outside the system con- 
sidered, the effects of cosmological sources. — For linear 
vorticity perturbations of FRW cosmologies no boundary 
conditions of the encoding type are needed because of the 
exponential suppression of super- H-dot scales. This re- 
sult is in contrast to Einstein's conclusion that the prob- 
lem of inertia and of boundary conditions could only be 
solved if the universe is spatially closed. 

In sect. XII we add scalar and tensor perturbations, 
and give the mathematical expression for Mach's princi- 
ple in this more general case. — We discuss Einstein's for- 
mulation of Mach's principle in 1918 |9j. For linear cos- 
mological perturbations Einstein's formulation is (1) too 
strong to be valid, as has been stated many times (grav- 
itational waves can exist without matter), and (2) un- 
necessarily strong, as we shall explain, for Mach's orig- 
inal purpose (the dragging of axis directions of inertial 
frames) . 

In sect. XIII we list open problems in connection with 
Mach's principle. 

In a companion paper p3| we shall present the results, 
analogous to those in this paper, for perturbations of 
FRW universes with k — ±1, i.e. spatially spherical or 
hyperbolic. — We use the conventions of Misner, Thorne, 
and Wheeler 0]. 

III. VORTICITY PERTURBATIONS 

Linear cosmological perturbations, see J. Bardeen [2l| . 
Kodama and Sazaki [25j], decouple into three sectors, 3- 
scalars (density perturbations), 3- vectors (vorticity per- 
turbations), and 3-tensors (gravitational waves). In the 
vector sector all quantities must be constructed from 3- 
vector fields with vanishing divergence. 

There are two crucial results, which make the analysis 
of vector perturbations very simple: 

1. There is no ambiguity about the slicing of space- 
time in the vector sector: Since the 3-scalar Sgoo 
must be zero in the vector sector, the lapse function 
is unperturbed. This was already emphasized by 
Bardeen |2l|: "The gauge-transformation proper- 
ties of vector perturbations are much simpler than 
for scalar perturbations, since now there is no gauge 
ambiguity about the time coordinate." 

2. The intrinsic geometry of each E t remains unper- 
turbed for linear vector perturbations. To prove 
this, one starts from the general metric in the vector 
sector, and one makes a gauge transformation pf| 
such that the 3-metric is unperturbed. — Hence for 
vector perturbations of FRW with k = {0, +1, — 1} 
the fixed-time slices E t remain Euclidean 3-spaces 
E 3 , resp 3-spheres S 3 , resp hyperbolic 3-spaces H 3 . 

Since linear vector perturbations leave 3-space in FRW 
unperturbed, the natural choices of spatial coordinates 
are comoving Cartesian coordinates for Euclidean 3-space 



in FRW with k = 0, resp comoving spherical coordinates 
for 3-space in FRW with k = {0, ±1}. The line elements 
are 

ds 2 = -dt 2 + Yy{ahi) 2 {dx 1 ) 2 + 2f3 l dx i dt], (4) 

i 

where (3 l is the shift 3-vector, the only quantity referring 
to vorticity perturbations in our coordinate choice, and 
(3i = {ahi) 2 f3 l . For a FRW background with k = and 
for E t in comoving Cartesian coordinates 

hi = l, 2 = 1,2,3. (5) 

For a FRW background with k = {0, +1, — 1} and for E t 
in comoving spherical coordinates {x, 0, 4>} 

h x = l, h g = R(x), h$ = R(\) sin 6, 

R{x) = {x> sinx, sinhx}- (6) 

In sections III, IV, V we give unified computations for 
FRW with k = {0, ±1}, but afterwards in this paper we 
specialize to k = 0. 

For perturbed FRW with k = all E t retain a Eu- 
clidean 3-geometry. With our choice of comoving Carte- 
sian coordinates for Euclidean 3-space the 3-coordinatc 
lines are geodesies in 3-space. In all other gauges, where 
a perturbed 3-metric g^ is used for Euclidean 3-space 

(in Bardeen's notation with 0), the 3-coordinate 

lines will quickly wind up more and more into spirals 
relative to geodesies on E t . Hence all other gauges (par- 
ticularly time-orthogonal) give a very awkward way to 
coordinatize Euclidean 3-space. 

We consider vector perturbations in an asymptotically 
unperturbed FRW universe, i.e. (3 l — > for r — > oo. We 
fix our coordinates to quasars in the asymptotic unper- 
turbed FRW space. (This is analogous to "star-fixed co- 
ordinates" for the physics of rotating black holes.) This 
fixes our gauge completely. — The basis vectors in the co- 
ordinate basis, <3i(P) = d/dx z , point along geodesies in 
E t from P to fixed quasars in the asymptotic FRW space. 
The case of perturbed FRW with k = is particularly 
simple, because 3— space remains Euclidean everywhere, 
hence each E t possesses global parallelism. 

It will turn out in Eq. 127|) that our shift vector j3 is 
equal to the gravito-magnetic vector potential A g , 

(3 = A e . (7) 

Our sign convention for the shift vector (3 Z is dictated 
by our wish to have no minus sign in Eq. Q). Our sign 
convention for the shift agrees with Misner, Thorne, and 
Wheeler (24|, and it is the opposite of Bardeen's |2J. — 
With our gauge-fixing the gravito-magnetic vector po- 
tential A g is directly related to a measurable geometric 
quantity: div A g — 0, and curl A g = — 2 r2 gyro , where 
r^gyro is the precession of local gyroscopes relative to the 
geodesies to asymptotic quasars. See Eqs. (|11H and l|28|) 
below. 
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Bardeen's gauge-invariant combination "J is defined by 
# = B {1) - \H^ ] (in Bardeen's notation). is the 

shift amplitude, and is the amplitude of the per- 

turbation of the 3-metric chosen for the unperturbed 3- 
geometry. Evidently the gauge-invariant amplitude ^ is 
equal to the shift amplitude B^ in the gauge which is 
fixed by = 0, i.e. in the gauge, where ^gij is un- 
perturbed. Bardeen emphasized in [26(: "Much has been 
made by some authors of the advantages of using gauge- 
invariant variables. However ..." We now transcribe 
Bardeen's formulation j2^| from the scalar sector to the 
vector sector: "The advantages of the gauge-invariant 
variable '5 are the advantages of working with the shift 
amplitude BW in the gauge where ^gij is unperturbed, 
no more and no less." 

The fact stated in the preceeding sentence expressed 
in our notation: The shift vector (3 l in the Cartesian 
gauge (for vector perturbations of FRW with k = 0) is 
equal to Bardeen's gauge-invariant variable 'J times his 
vector harmonic Q^ l (x) (for a given wave number and 
polarization) apart from the sign discussed above, 

p i (x,t)=Ai = -V(t)QW i (x). (8) 

Our gravito-magnetic vector potential A g is directly pro- 
portional to Bardeen's gauge-invariant combination "P. 
Working with our gauge, i.e. with comoving Cartesian 
(or spherical) coordinates for Euclidean 3-space, is totally 
equivalent to working with Bardeen's gauge-invariant po- 
tential ty. 



IV. FIDUCIAL OBSERVERS AND 
GRAVITOMAGNETIC FIELD 

Our aim is to obtain the laws of linearized gravito- 
magnetism in a formulation analogous to electromag- 
netism with a split to (3 + l)-dimensions. What is the 
operational definition for the gravitoelectric field E g and 
gravitomagnetic field B g ? According to the equivalence 
principle for a free-falling, non-rotating observer there are 
no gravitational forces at his position, E g = 0, B g = 0. 
It all depends on the choice of fiducial observers (FIDOs) 
with their local ortho-normal bases (LONBs). see Thome 
et al |27j. For FIDOs which are not free- falling and/or 
which are rotating (relative to local gyroscope axes) there 
are gravitational forces (inertial forces, fictitious forces, 
e.g. the Coriolis force). 

The operational definitions of E g and B g are inde- 
pendent of perturbation theory. They involve a fam- 
ily of world lines of fiducial observers, FIDOs, of any 
given choice, such that through each space-time point 
P there passes precisely one world line, i.e. the fam- 
ily of world lines forms a congruence. The measure- 
ments of a FIDO are directly encoded by the compo- 
nents (of momentum p and of gyroscope spin S) with 
respect to the FIDO's local orthonormal basis, where 



eg(P) = ufido(-P), ancl ej(P) are given by any definite 
choice. Hats always refer to LONBs, and bars designate 
space-time vectors. — To convert local measurements by 
a FIDO to LONB-components p a and S a one only needs 
r; a g = diag{ — 1, +1, +1, +1}, but one does not need <? M „, 
which is not available before having solved Einstein's 
equations. Therefore it is necessary to abandon vector 
components in a coordinate basis, to work with LONBs, 
and to use the formalism of Elie Cartan j^, |2^. — 
Measuring {E g , B g } involves FIDOs measuring first time- 
derivatives along their world lines, on the one hand of the 
momentum components pi of free-falling quasistatic test 
particles, and on the other hand of the spin components 
S- of gyroscopes carried along by the FIDOs, 

— p\ = mEf free-falling quasistatic test particle, (9) 
^5- = - i[J5 g x S\i gyro comoving with FIDO, (10) 

where t is the local time measured by the FIDO. Arrows 
denote 3-vectors in the tangent spaces spanned by the 
spatial legs of LONBs. E g = g is the gravitational accel- 
eration of free-falling quasistatic test particles relative to 
the FIDO. — The operational definitions of Eqs. (|9I1U|) 
are the same as for a classical charged spinning test par- 
ticle in an electromagnetic field except that the charge 
q is replaced by m, and the gyromagnetic ratio q/(2m) 
is replaced by 1/2. Eq. (|10fl gives the angular velocity 
of precession of the gyroscope's spin axis relative to the 
axes of the FIDO, 

^gyro = — 2^6' (11) 

which is an equivalent operational definition of B g . — For 
gyroscopes carried along by FIDOs the magnitude of S 
is constant, and for non-rotating FIDOs the components 
of the vector S are constant, because a gyroscope ex- 
periences no mechanical torques by construction and no 
gravitational torques by the equivalence principle. 

From the operational definitions of {E g ,B g } it follows 
immediately that these fields arc identical (apart from 
an overall sign) to the connection coefficients (w a j,)g for 
a displacement along the FIDO's world line. — The con- 
nection 1-forms Co a ~ resp. their components in LONBs (= 

Ricci rotation coefficients), (ui a t)c, are defined by 

V a e g = e a (w e s ) a . (12) 

In words: the Ricci rotation coefficients (w c g)a give the 
rotation resp. the Lorentz boost (^ c g) of the LONBs 
under parallel transport along eg. Parallel transport of 
u = eg along u is given by free fall, parallel trans- 
port for is given by gyroscope axes. Relative to FI- 
DOs the equations of motion for free-falling test parti- 
cles (geodesic equation) and for spin axes of gyroscopes 
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(Fermi transport) specialized to gyroscopes carried along 
by a FIDO are 

^ + K)s ^ =0 , f + (^=0. (13) 

With Eqs. lfT3|) the operational definitions Eqs. i|9llTTl) are 

translated into the equivalent definitions involving con- 
nection coefficients with a displacement index 0, namely 
a Lorentz boost <x>jg per unit time (acceleration g = E g ) 
resp a rotation angle per unit time (angular velocity 

(ux) 6 = -Ef, (^o^—Bfj, (14) 

where /;, , = e^-Bj, fiy = e«j%, and <>,, = (w-.)o- Thc 
last equation (angular velocity) gives the motivation for 
using the letter u> for Ricci rotation coefficients (connec- 
tion coefficients). Only displacments along eg appear in 
Eq. (| 1 4p . because for quasistatic test particles (typically 
initially at rest at the position of the FIDO) the spatial 
separation from the FIDO grows only quadratically in 
time. 

Our choice of FIDOs: The world lines of our FIDOs 
are at fixed x 1 in our coordinates, which are fixed to 
quasars in the asymptotic FRW universe. The 3- velocity 
of our FIDOs relative to the normals on E t is equal 
to the shift 3- vector (3 % in Eq. (0J. — The construction 
of the world-lines of FIDOs is equivalent to the gauge- 
invariant statement that our FIDOs are at rest in the 
Hubble frame given by asymptotic quasars on E t . — We 
choose the spatial basis vectors of our FIDOs, ej(P), fixed 
to directions of geodesies on E t from P to quasars in the 
asymptotic FRW universe, again a gauge-invariant state- 
ment. Hence we fix ej(P) in the same spatial directions 
as e~i(P) = d/dx l , i.e. in 4-space the directions of and 
e, differ by a pure Lorentz boost, see Eq. lfT3|l below. 

V. CONNECTION COEFFICIENTS AND 
EQUATIONS OF MOTION FOR MATTER IN 
GRAVITOMAGNETISM 

Cartan's formalism with local orthonormal tetrads is 
crucial for gravitomagnetism, because (1) the operational 
definitions of the gravito-electric and gravito-magnetic 
fields (independent of perturbation theory) are given by 
measurements of the corresponding FIDOs, Eqs. © - 
(|llfl . and (2) the measured LONB-components of the 
gravito-electric and gravito-magnetic fields are identical 
to connection coefficients in the local orthonormal ba- 
sis of the FIDO, Eqs. l(Ti|> . A third reason, why it is 
necessary to work in Cartan's formalism with local or- 
thonormal tetrads, will be presented in section IX. 

The computation of connection coefficients in Cartan's 
formalism: We give a short, self-contained derivation 
of the formulae needed from Cartan's formalism with 
LONB-tetrads. This should be helpful for readers only 



familiar with the literature on cosmological perturbation 
theory and with standard one-year graduate courses on 
General Relativity. Readers not interested in calcula- 
tional methods can go directly to the results starting with 
Eq. (23. 

A. Cartan's first equation in LONB components 

The first step is to express our choice of LONBs e a (P), 
given at the end of the previous section, in terms of the 
coordinate bases e a (P) = d/dx a , i.e. e"a = (ea) Q e a . To 
first order in (f3 l /c) 

eg = e , e k = (efe + /9fee ). (15) 

Since the spatial LONBs point in the same spatial direc- 
tions as the spatial coordinate bases, we use latin letters 
from the middle of the alphabet both for spatial LONBs 
(with hat) and for spatial coordinate bases (without hat). 
The dual bases (basis 1-forms) 9 a for LONBs are de- 
fined by (9 a ,ef ) ) — <5?, where tildes designate space-time 

1-forms. The LONB 1-forms 9 a are expanded in the co- 
ordinate basis 1-forms, 9 a = dx a , i.e. 9 a = (9 & ) a 9 a , 
by 

9 6 = 9° - (3 k 9 k , 6 % = ah k 9 k . (16) 

The coefficients of the inverse expansion, i.e. coordi- 
nate bases in terms of LONBs, are (e a ) a = (9 a ) a resp. 

(9 a )a = (ea) a . 

The second step is computing the exterior derivative d 
of the basis- 1-forms, {d9 c ) a j3 = d a (9 c )p — dp(9 c ) a , where 
[a/3] must be in the coordinate basis. Then one converts 
to components [ab] in the LONB, 

(d9") &h ^ -C.J = (e & r[d a (e% - d p (9%](e b f. (17) 
The coefficients C - c are identical to the commutation 

ab 

coefficients of the basis vectors, 

[e a ,e s ]^C aS £ e e . (18) 

This is easily shown by using d a {6 c ,e<£) — = 
[d a (9 & ) f ,](e- b r + (9 & ) f3 [d a (e i n 

The third step is obtaining the connection coefficients 
from the commutation coefficients. The definition of the 
connection via basis 1-forms is (y a 9 c )p — ~{uo c ^) a {9 d )(j. 
We take both the displacement index a and the equa- 
tion's component index (3 in the coordinate basis, and 
we antisymmetrize in [a(3\. This makes the Christoffel 
symbols T p a/3 on the left-hand side disappear, since they 
are symmetric in a, (3. Hence the left-hand side reduces 
to (d9 c ) a p. Dropping the equation's component indices 
[a, (3] gives 

d¥ = -Q & < A ¥. (19) 
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This is Carton's first equation. The wedge product (ex- 
terior product) of two 1-forms a and p is (a A p) a p = 
o~ a P0 — o-pp a . Taking Cartan's first equation in LONB 
components gives the commutation coefficients, Eq. (|17|l . 
on the left-hand side, and the right-hand side simplifies in 
LONB because {d d )~ b = 6$. Hence Cartan's first equation 
in LONB components is 

c d = K)« - (<4)s- (so) 

This equation is easily solved for the rotation coefficients, 

(LU.,)a = -fC-f- +C.~*-C*~A. (21) 
v cb> a 2 c " a ca " " ac 



in the tangent spaces to S t . From Eqs. I|23l I24|) . we see 
that the shift vector f3 must be identified with the grav- 
itomagnetic vector potential A g . From Eqs. 1231 12^|) 
follow 

g oi = fa = (A g ), (27) 
B g = curL4 g , E & = --d t (aA g ), (28) 

curing + -^d t {a 2 B g ) = 0. (29) 

Equations H28fl and l|29|) are identical with the homoge- 
neous equations for electromagetism in FRW space-times 
with k = 0,±1. 



B. Connection coefficients for vorticity 
perturbations on a Minkowski background 

To first order in j3 l and with Cartesian spatial coordi- 
nates on E t , the commutation coefficients are very simple 
to compute, because only (8°)i = —f3i is space-time de- 
pendent, and the prefactors in Eq. (|TT|) . (ea) Q and (e^)' 3 , 

can be set to 1. Hence C.~° = (d(3) a b, and 
(^lo)o = - Ef = dtpi, 

= -\b^ = -\ma = («a)j. 

(u%)t = 0. " (22) 

All connection coefficients with respect to LONB's are 
directly measurable (in contrast to Christoffel symbols, 
which refer to coordinate bases). 

C. Connection coefficients for vorticity 
perturbations on FRW with {k = 0, ±1} 

It is again straightforward to compute the commu- 
tation coefficients and the connection coefficients using 

Eqs. iini csi czi I2U) , 

(a, ;e ) 6 ee -Ef = -d t {a^), (23) 
1 a 

N)o - -\Bl = -\w)r 3 , (24) 

(«m) 3 = k; • •>".,"• (25) 

(f»13h = hk( H l 3 3 + d 3 L i)- S 'd H 0i + d i L ' j )' ( 26 ) 

where Lj = log/ij, <9; = (ah i )~ 1 di, and i? is the Hubble 
rate. In Eq. (|26|l there is no summation over repeated 
indices on the right-hand side. — The Hubble term in 
Eq. I|25|) follows directly by comparing the operational 
definitions of the Hubble rate and of the connection co- 
efficients. 

Since we work to first order in the vorticity perturba- 
tions (i.e. in /3 l ), we can identify E g and B g with vectors 



D. Equation of motion for free- falling test particles 

The equation of motion (geodesic equation) for free- 
falling test particles of arbitrary velocities v < c (e.g. 
photons) in linear vorticity perturbations on a Minkowski 
background reads 

±( Pi )=e[E g + (vxB g )} l . (30) 

This is identical with the Lorentz law for electromag- 
netism, except that the charge q is replaced by the en- 
ergy e of the test particle, e = mo(l — ft 2 )" 1 / 2 for massive 
particles, and e = hv for photons. With Eq. and in 
a stationary gravitomagnetic field (E g = 0), Eq. i|3U|) 
becomes 4x{p%) = — 2e[u A f2 gyro ]-, the Coriolis force 
law. Note that £l gym is minus the rotation velocity of 
the FIDO relative to the gyroscopes' axes. A homo- 
geneous gravitomagnetic field can be transformed away 
completely by going to FIDOs which are nonrotating rel- 
ative to the spin axes of local gyroscopes, i.e. in homoge- 
neous rotation relative to the previous FIDOs. Physics 
in a homogeneous gravitomagnetic field is equivalent to 
physics on a merry-go-round in Minkowski space. — The 
centrifugal force is quadratic in fi, therefore it is missing 
in linear perturbation theory. 

For vorticity perturbations on a FRW background with 
k = in Cartesian comoving LONBs we obtain 

^j t (a Pi )=e[E g + (vxB g ) + Hvx (|x% (31) 

For FRW with k = {0, ±1} in spherical comoving LONBs 
there are additional terms from c^Lj in Eq. (|26|l . These 
terms are present even in the absence of vorticity pertur- 
bations and of Hubble expansion, because in a spherical 
basis for k = spatial LONBs are not parallelized, and 
for k = ±1 global parallelization is impossible on St. 

VI. EINSTEIN'S EQUATIONS FOR 
GRAVITOMAGNETISM: AMPERE'S LAW 

The computation of curvature in Cartan's formalism: 
We give a short, self-contained derivation of the formu- 
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lae needed from Cartan's formalism. Readers not inter- 
ested in calculational methods can skip the first part of 
this section and go directly to the results starting with 
Eq. JS§. 

A. Cartan's second equation in LONB components 

The curvature 2-form 1Z a ^ has LONB components 

(7\l a g) £ j, which are the LONB components of the Rie- 

mann tensor R a ~ «. The Riemann tensor can be opera- 
oca 

tionally defined by the action of (V 7 Vi — ViV 7 ) on the 
LONB 1-forrn 9 & , 

(V 7 V 5 - V«V 7 )£ a = -9 b (K% s , (32) 

where the covariant derivatives V 7 and V5 must be in the 
coordinate basis. To compute the curvature 2-form we 
first use Vs9 a = — (uj a .)s6 e . Acting on the right-hand side 
by V 7 gives two terms. One term comes from V 7 acting 
on 9 e , and after antisymmetrization in [jS] it produces 
— (uj a g A u e ^) 1 $9 b . The other term comes from V 7 acting 

on the expansion coefficient (number field) (u> a ~)s, where 
it can be replaced by c? 7 , and after antisymmetrization 

in [y5] it produces — (du! a ^)~ f $9 b . Hence we obtain 

1l% = duj\ + u>%Acj% (33) 

which is Cartan's second equation. 

To obtain the LONB components of the first term of 
the right-hand side, (dco a ^) £ ^, we must first convert the 
connection components of Eqs. (|23I - I26|) from the LONB 
to the coordinate basis, (uj a ^)s = {uj a ^ c i{9 d ) s , then take 

the exterior derivative, d~,{(uj a -)j(0 d )s} — [7 <-> 8\, and 
then convert the result back from coordinate compo- 
nents to LONB components. The partial derivative of 
the product gives two terms, one with a partial deriva- 
tive of (w a g)j, the other with d 1 {9 d )$, which produces 
another connection- 1-form component. In the second 
term of Cartan's second equation these conversions from 
LONB to coordinate basis and back again cancel, since 
there is no derivative in between. The result is Cartan's 
2nd equation in LONB components, 

= Keara 7 (4) d --(4) ; (a;4)e + (<)a(4) d -] 
-[c<->4 (34) 

B. Einstein equations for vorticity perturbations of 

Minkowski space 

For linear vorticity perturbations of Minkowski space 
(with Cartesian coordinates for 3-spacc) all non-zero con- 
nection coefficients in Eqs. (1231 - 12611 are of first order in 
the perturbations. Therefore the second term of Car- 
tan's second equation can be neglected, and in the first 



term one need not distinguish components in LONB from 
components in the coordinate basis. For vorticity pertur- 
bations the important Einstein equation is the equation 
for Gg; = R bh 

% = = (^oj)lj = (35) 

Hence Einstein's Ggj-equation for general, i.e. time- 
dependent, vorticity perturbations on a Minkowski back- 
ground reads 

curLB g = -16nG N J E . (36) 

In contrast to the Ampere-Maxwell equation of ordi- 
nary electrodynamics, the Maxwell term (dtE) is ab- 
sent in the time- dependent, i.e. non-stationary context 
of gravitomagnetodynamics. Eq. I)36|) is consistent with- 
out (dtE g ), because in the vorticity sector all vector 
fields are divergence-free, particularly divJ e = 0. A term 
dEg/dt cannot be present in gravitomagnetism, other- 
wise Eq. (|36l) together with Faraday's law for gravito- 
magnetism, Eq. H29|) . would erroneously predict gravi- 
tational vector waves. — Eq. I|36|) for gravito-magneto- 
dynamics (i.e. for a non-stationary context) is identical 
to the original law of Ampere for stationary magnetism, 
except that the charge current J q is replaced by the en- 
ergy current J e , and the prefactor 47r is replaced by the 
prefactor (— I67tGn)- — The source J\ = T° l is the en- 
ergy current density, which is equal to the momentum 
density (for c = 1). In linear perturbation theory for a 
perfect fluid J e = (p + p)v. 

The Gjq equation is an equation at fixed time, a con- 
straint equation, called momentum constraint, since the 
momentum density appears on the right-hand side of 
Eq. To see the analogous structures of gravit- 

omagnetism and electromagnetism, it is more instruc- 
tive to formulate this constraint equation, as we have 
done in Eq. H3(jfl . y i a connection 1-forms, which in- 
volves (difij — djf3i), i.e. the gravitomagnetic field, than 
via the extrinsic curvature tensor Kij, which involves 
(di(3j+djPi). Of course the resulting constraint, if written 
in terms of A g — j3, is the same, AA g — 16ttG^J £ . 

The <5Ggg equation with the source STq^ is trivially 
fulfilled, since these two objects are 3-scalars and there- 
fore vanish in the vector sector. — The source 5T:~- van- 
ishes (for linear vorticity perturbations and perfect flu- 
ids), since it is of second order in the perturbation. The 
(5G-- equations give do(dif3j +dj(3i) — 0, i.e. the shear of 

the field f3 has vanishing time-derivative. 



C. The momentum constraint for vorticity 
perturbations of spatially flat FRW spaces 

With Cartesian comoving coordinates for flat 3-space 
there are two new terms in the connection coefficients, 
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(^o)f W = h'o H and K)j = - Com- 

puting = (^oj)ij Eq- P^! 1 we obtain the corre- 
sponding Einstein equation, 

curLB g - 4{dH/dt)A g = -16ttG n J s , (37) 

where we have used j3 = A s . The scale factor a of the 
spatially flat FRW universe does not appear in this fixed- 
time equation. Note that 'fixed-time equation' means 
that there are no partial time-derivatives of the fields 
B g {x,t), etc. The coefficient (dH/dt) is merely a given 
input number at each time. — Every symbol in this equa- 
tion refers to the physical scale. This is a general fact: 
If we consider any fixed-time equation at any time t±, we 
are free to set a(ti) = 1, hence the scale factor disap- 
pears. 

Friedmann's equations give ff-dot, 

dH/dt = -4ttG n ( P + p). (38) 

Since {dH/dt) < for {p/p) > -1, we define the ff-dot 
radius RdH/dt by R^n/dt = {—dH/dt) -1 , and we define 
M by 

p 2 = -4{dH/dt) = 4{R dH/dt )- 2 . (39) 

We use B g — curlA g , and with divA g = we have 
(curl curl A g ) = —AA g . Therefore Eq. (|37|) becomes 

(-A + p 2 ) A g = -16 7T G N J e , (40) 

an elliptic equation. Concerning the absence of par- 
tial time-derivatives of A g in Eq. (|40|l for the time- 
dependent context of gravitodynamics, see the comments 
after Eq. (|36[) . — The new term on the left-hand side, 
{—AH A g ) = {p 2 A g ), is responsible for the exponential 
cutoff on scales larger than the iJ-dot radius, Eqs. QEJ, 
which is crucial to remove the problem of 'overdragging ', 
as shown in Fig. 1, and for obtaining a weight function 
with normalization to unity in our Eq. which gives 

figyro as the weighted average of ^matter- 

Bardeen |2l| gives two alternative choices for gauge- 
invariant amplitudes for the matter velocity, v s and v c 
in his Eqs. (3.21, 3.23). — Bardeen uses Einstein's mo- 
mentum constraint in coordinate- basis components, Tj 
on the source side, and this is directly proportional to 
v c . — In contrast we use the momentum constraint in 
LONB components, Tg-., and this is directly proportional 
to Bardeen 's other gauge-invariant variable v s . 

The physical meaning of v s and v c : The field with 
the amplitude v c , used by Bardeen, is directly related 
(via its curl) to the angular velocity field of matter rela- 
tive to axes of local gyroscopes all over H t . But this has 
the severe drawback that v c is not measurable without 
prior knowledge of the solution of Einstein's equations, 
g^u and therefore A g and B g . Hence v c , employed by 
Bardeen, cannot be used as a matter input for solving Ein- 
stein's equations in the context of Mach's principle and 



of dragging, where a measured matter input is needed 
before having solved Einstein's equations. See also sects. 
VIII and IX. — In contrast we use Bardeen's other gauge- 
invariant amplitude v s , which is directly related to the 
angular velocity of matter relative to the asymptotic un- 
perturbed quasars. This is an input which is measurable 
without prior knowledge of got = A^ g \ i.e. before knowing 
the output of solving Einstein's constraint equation. 

The difference between the two gauge-invariant ampli- 
tudes for the matter velocity is given by (v c — v s ) = — ty, 
Eq. (3.23) of Bardeen. The ^-term, a geometric term, 
which is on the right-hand side of Bardeen's momentum 
constraint, must be moved to the left-hand side, the geo- 
metric side. With the prefactors in the momentum con- 
straint and in our notation the 'I'-term appears as p 2 A g . 

Both sides of our Eq. (|40|l are gauge-invariant, the 
left-hand side because our A g is directly proportional 
to Bardeen's gauge-invariant potential Eq. (jJJ, the 
right-hand side because our matter velocity v is directly 
proportional to Bardeen's gauge invariant amplitude v s . 

VII. MACH'S PRINCIPLE 

A. Evolution of inertial axes determined 
exclusively by cosmic matter flows 

The solution of Eq. I|40l) is the Yukawa potential for 
A g = (3 in terms of the energy currents J £ at the same 
time, 

A g {r,t) = -4G N / d 3 r'J e (?,t) g^L^Z r 3 . (4 i) 
J \r — r'\ 

This is analogous to the formula for ordinary magne- 
tostatics except for the exponential cutoff, i.e. the 
1/r potential is replaced by the Yukawa potential 
(l/r)exp(— /xr). The Green function which is exponen- 
tially growing for r' — > oo has been rejected on the stan- 
dard grounds of field theory. This gives our first conclu- 
sion: The contributions of vorticity perturbations beyond 
the _ff-dot radius are exponentially suppressed. 

The fundamental law of gravitomagnetism in integral 
form for linear vorticity perturbations and for spatially 
flat FRW universes {k = 0), the equation for B g and 
rig yro in terms of sources at the same time, is 

B g {P) = -4G N j d?r Q [n PQ x J e {Q)] Y^{r PQ ) (42) 
-d e~^ r 

Y„ (r) = — = Yukawa force, (43) 

dr r 

n gym = -is g (r gyro ), (44) 

where npQ is the unit vector pointing along the geodesic 
on E t from P to Q, and Q stands for "source point" . We 
have applied B g = curlA g to Eq. I|41|) . 
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Both B g = — 2f2 gyro and the transverse velocity in 
[ripQ x J e (Q)] are measured relative to geodesies on E t 
from the the gyroscope at P to quasars in the asymp- 
totic FRW universe. Since £t is a Euclidean 3-space, 
global parallelism on S t can be used when comparing 
vectors at P and Q. Although we use Cartesian coor- 
dinates for Euclidean 3-space, all quantities in Eq. I|42|) 
have a gauge invariant meaning, they are directly mea- 
surable, e.g. d 3 rg is the measured volume element, and 
rpQ is the geodesic distance from P to Q. — The scale 
factor a of the spatially flat FRW universe does not ap- 
pear in the fixed-time equation l|4*2")> . as discussed after 
Eq. l|57|) : every symbol refers to the physical scale. 

The only differences from Ampere's law in integral 
form are: 

1 . the replacement of the current of electric charge J q 
by the measured energy current of matter J e , 

2. the factor (— Gn), as in the transition from 
Coulomb's law to Newton's law, 

3. the additional factor 4, which occurs in the transi- 
tion from Ampere's law of ordinary magnetism to 
gravitomagnetism, 

4. the replacement of the 1/r 2 force in Ampere's law 
by the Yukawa force with its exponential cutoff, 
Eq. H43|) . which occurs in the transition from a 
Minkowski background to gravitomagnetism on a 
a FRW background with k = 0. 

An analogous Yukawa cutoff at super- H-dot scales in 
the scalar sector of cosmological perturbation theory oc- 
curs in the 'uniform expansion gauge', which is also called 
the 'uniform Hubble gauge'; see Ref. |29|. 

Note the fundamental difference between cosmological 
gravitomagnetism, Eq. I|42l) . and Ampere's law. The lat- 
ter does not hold in a rotating reference frame, unless 
one introduces fictitious forces. The same is true for the 
equations for electromagnetism in special relativity, and 
for general relativity of the solar system in asymptotic 
Minkowski space, where fictitious forces must be encoded 
by boundary conditions at spatial infinity as explained in 
Sec. XI. In contrast, for cosmological gravitomagnetism, 
Eq. I|42l) remains valid, as it stands, in a frame which is 
rotating with angular velocity f2* relative to asymptotic 
quasars, as explained in sec. II after Eq. This es- 
tablishes the fact that the asymptotic inertial frame has 
no influence in cosmological gravitomagnetism. The lo- 
cal nonrotating frame (at any point P) is not needed as 
an input in applying Eq. (|42|l . i.e. no absolute element is 
needed as an input. The time evolution of local inertial 
axes is an output, determined exclusively by the weighted 
average of cosmological energy flows. 

Eqs. f4*2")) - l |4"4*jl state what specific average of the mea- 
sured energy flow of matter out there in the universe 
determines the motion of gyroscope axes here. This an- 
swers Mach's question |30j : 



"What share has every mass in the deter- 
mination of direction ... in the law of inertia? 
No definite answer can be given by our expe- 
riences." 

The gravitomagnetic moment density p g per volume 
element is the analog of the magnetic moment den- 
sity (with J q replaced by the measured energy current 
J £ ). The gravitomagnetic moment density appears as 
the source in Eq. (|4*2|) . and it is equal to one half 
of the measured angular momentum density of matter 
L = (p + p) (f x v) at the source point Q relative to the 
gyroscope point P, 

i«g = \[rpq X J e {Q)] = \l pq , (45) 

^gyro(P) = ~\b & {P) 

= 2G N f d 3 r Q ^ Y^(r PQ ). (46) 
J r PQ 

It is a matter of taste and choice of emphasis, whether 
the term 'source' for the gravitomagnetic field is used for 
(1) the measured energy current (= momentum density) 
J e as in the analogue of Ampere's law (momentum con- 
straint), Eqs. H36|) and JSJJ, or (2) the measured angular 
momentum density L = (r x J £ ), as in the integrated 
form of the analogue of Ampere's law, Eq. I|46(l . 

B. Evolution of inertial axes exactly follows the 
average of cosmic matter flows 

The gravitomagnetic moment density per r-interval, 
dp g /dr, is the lowest multipole term, the £ = 1 term, of 
the "odd parity sequence" , P = (— 1) +1 , in the multipole 
expansion of the vector source J e for r ODS = ''gyro = and 
''source > 0. Other multipoles cannot contribute to the 
gyroscope's precession, i.e. to the gravitomagnetic field 
B g at r = 0, for reasons of symmetry under rotations 
and space reflection. — In first-order perturbation theory 
(p + p) ^constant all over space, hence also the velocity 
field at a given radius r can only contribute via the term 
with {£ — 1, odd parity sequence}, which is equivalent 
to a rigid rotation with the angular velocity ^matter (?*)■ 
Using this fact and Eqs. 

we obtain 

^gyro = — / dr r O mat ter(r-) G M (r), (47) 
Jo 

<7 M (r) = r 2 Y^r) = (1 + pr) cxp(-pr). (48) 

This is our most important result. Eq. 1|47JI shows that 
f2gy ro is the weighted average of ^matter, i-e. the preces- 
sion of a gyroscope, ^ gyr o, exactly follows, i.e. it is exactly 
dragged by the weighted average of ri mat t e r- — The weight 
function is normalized to unity, as it must be for an av- 
eraging weight function in any problem. This crucially 
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depends on the exponential cutoff in the Yukawa force 
Y^(r). The cutoff function GAr) is 1 for r -C and 
goes to zero exponentially fast for r S> jU — 1 . The weight 
function per logarithmic r- interval in Eq. (|47|l is r 2 C /i (r), 
i.e. it grows quadratically until one reaches, roughly, the 
i?-dot radius, and then it goes to zero exponentially. — 
The conclusion is that Eq. I|47|) is a clear demonstration 
of how Mach 's principle is contained within cosmological 
general relativity. 

When talking about the "angular velocity of cosmic 
matter" one must keep the following in mind: For a 
given arbitrary velocity field the densities of gravitomag- 
netic moment and of measured angular momentum at 
one source point Q relative to the gyroscope point P 
are well-defined. But one cannot uniquely define a den- 
sity of angular velocity at one source point Q relative to 
P, P), because one can add terms r P Qf(Q) without 

changing v(Q). 

C. The dragging fraction 

The dragging fraction refers to the simple cosmologi- 
cal model discussed in the Introduction, a uniform (rigid) 
rotational motion with Sl ma tter constant inside a radius 
i? ro t around a gyroscope and zero outside and with en- 
ergy density and pressure spatially constant everywhere. 
The definition of the dragging fraction /drag is given by 

^gyro — /drag ^matter ? (^^) 

where ri gyr o and f2 ma tter are measured relative to 
geodesies on E t from the gyroscope to quasars in the 
asymptotic FRW space. — A better word than dragging 
fraction is voting power of the rotating matter (in con- 
trast to the voting power of the non-rotating matter 
outside), because there is exact dragging of ^ gy ro by 
a weighted average of ri ma ttor according to Eq. I|47(l . — 
Putting fi m attcr = constant inside a radius R TOt and 
^matter = outside in Eq. (|4T|) we obtain the dragging 
fraction, 

p 2 f Rrot 

/drag(Prot) ^yj dr r (1 + fir) exp(-pr) 

= 1 - {cxp(-H [1 + p,r + fJ?r 2 /3]} r = Rtot . (50) 

The dragging fraction (voting power of the matter inside 
i? ro t) is shown by the solid curve in Fig. (1). It grows 
quadratically with i? ro t until it reaches order of one for 
-Rrot = RdH/dt- As R rot grows beyond R d H/dt, the drag- 
ging fraction approaches the exact value 1 exponentially 
fast, i.e. exact dragging of the gyroscope axes by the 
rotating matter. This holds for any equation of state. 

The problem of it 'overdragging', shown by the dashed 
curve in Fig. (1), arises when extrapolating the results 
of gravitomagnetic perturbation theory on a Minkowski 
background beyond the region of validity, as noted in 



[lij . — The solid line shows, how the problem of 'over- 
dragging' is removed by the exponential cutoff from 
cosmological perturbation theory for super-Hubble-dot 
scales, i.e. by the exp onential cutoff in Eqs. (|4UI 1421) . — 
Lynden-Bell et al |l6j do not have an exponential cutoff, 
because they use the conserved angular momentum in- 
stead of the measured angular momentum as the source. 
See our discussion below in sect. IVIIII 

Analogous anti-dragging effects in magneto statics: A 
rotating charged spherical shell acting on magnetic dipole 
moments inside has the opposite sign in Ampere's law 
compared to gravitomagnetism, Eq. (|36|l . and this causes 
antidragging in the magnetostatic case. 

Physics could have evolved differently after Mach 
(1883): Physicists at that time could have extended the 
correspondence between Coulomb's law and Newton's 
law and postulated gravitomagnetism in correspondence 
to Ampere's magnetostatics. Their prediction of the 
dragging of inertial frames would have been too small 
by a factor 4 compared to the correct result for linear 
gravitomagnetism on a Minkowski background. For the 
simple case of a homogeneous rotation out to some radius 
R of cold matter with p = constant they would have ob- 
tained perfect dragging for GM / (Re 2 ) — 2 and noted the 
problem of 'overdragging' for GM/ [Re 2 ) > 2. In the late 
1920's, for an approximately homogeneous and expand- 
ing universe with cold matter, this critical radius would 
have been identified with the Hubble radius (apart from 
a missing factor 2). 

Measuring everything relative to axes of gyroscopes at 
one given location (instead of relative to quasars in the 
asymptotic FRW space) makes the the left-hand side of 
Eq. H46|) vanish. Hence Eq. I|46|) reduces to the state- 
ment that the measured angular momentum of matter 
relative to gyroscopes at one given location will vanish 
after averaging with the weight r~ 3 and with the expo- 
nential cutoff C^(r). — From the point of view of mea- 
surements it is preferrable to mesure relative to gyro- 
scopes at one given location. But to sec the structure of 
gravitomagnetism and its correspondence with electro- 
magnetism most clearly, it is best to measure relative to 
quasars in the asymptotic FRW space. 



D. The de Sitter limit 

Dark energy with p/ p = —1, i.e. a cosmological con- 
stant, does not contribute in Mach's principle, Eq. I|42|l . 
since there is no flow of energy associated with it, its 
energy current J e — (p + p)v vanishes. — When consid- 
ering a FRW universe with k = and letting (p/p) get 
closer and closer to (—1) at fixed Hubble rate (hence 
fixed p) until we are arbitrarily close to a de Sitter uni- 
verse, Mach's principle continues to work (in linear per- 
turbation theory). When taking this limit the prefactor 
p? = (167tGn)(p+p) in Eq. I|5U|I gets smaller and smaller 
(and tends to zero), but the H— dot radius fx , which 
cuts off the otherwise quadratically divergent integral in 
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Eq. H50(l . gets larger and larger such that perfect drag- 
ging is maintained for arbitrarily small p: 2 , i.e. arbitrarily 
close to a de Sitter universe. 



E. Einstein's objection: Interaction of masses 

Einstein's objection to Mach's principle in his autobi- 
ographical notes of 1949 |3lT |: 

"Mach conjectures that inertia would have 
to depend upon the interaction of masses, 
precisely as was true for Newton's other 
forces, a conception which for a long time I 
considered as in principle the correct one. It 
presupposes implicitly, however, that the ba- 
sic theory should be of the general type of 
Newton's mechanics: masses and their inter- 
actions as the original concepts. The attempt 
at such a solution does not fit into a consis- 
tent field theory, as will be immediately rec- 
ognized." 

We have shown how this apparent difficulty is resolved 
in General Relativity, specifically in linear Gravitomag- 
netism: The relevant Einstein equation (the momentum 
constraint) has the form of Ampere's law with a Yukawa 
cutoff, Eq. H37fl . It is a remarkable fact that Einstein's lo- 
cal field equation for Gqj has the solution Eq. I|42|) . which 
has the form of an instantaneous action- at- a- distance. It 
says that the measured mass-energy flow out there in the 
universe does indeed determine the precession of gyro- 
scope axes here. — See also refs. 0] |22| . 



VIII. MEASURED MATTER INPUT FOR 
EINSTEIN'S EQUATIONS 

The input on the right-hand side of Eqs. (|42II47|> . which 
express Mach's principle, is the measured angular veloc- 
ity of matter (stars, galaxies, etc) relative to the geodesies 
on S t to quasars in the asymptotic FRW space. No 
knowledge of the metric perturbation (3 is needed when 
determining the input for Eqs. Q42I I47fl . The measured 
angular velocity of matter is a purely kinematic or ki- 
netic input. Similarly the measured angular momentum 
density L = (p + p) (r x v) = (p + p) (r 2 sin 2 6 fl) is 
a purely kinetic input. We call it the kinetic angular 
momentum. — The geometric-dynamical output of solving 
Einstein's equation is (3 = A g , the gravitomagnetic 
vector potential in Eqs. I|4(JI I41f> . hence the gravitomag- 
netic force B g in Eq. 1421) . — The measured, kinetic angu- 
lar momentum must be distinguished from the canonical 
angular momentum, which we introduce in the following 
paragraphs. 



Kinetic versus canonical momentum and angular 
momentum 

The Einstein- Hilbert action for linear vorticity pertur- 
bations (i.e. gravitomagnetism) on a Minkowski back- 
ground and for point particles gives 

S = (WttGn)- 1 J d^x^fg- (curll g ) 2 

+ / dt^[-mx 2 n + mx n A s (x n ,t)}. (51) 

It is natural to take the nonrelativistic expression [y <C c) 
for matter when discussing linear perturbations. Except 
for the prefactor — (167tGn) _1 in the first term and the 
prefactor m in the last term, the action is the same as 
for electromagnetism without the (dtA) 2 term. 

The Lagrangian gives the equations of motion via the 
standard Euler-Lagrange equations (as in classical me- 
chanics and classical electrodynamics) , if and only if the 
Lagrangian is defined by 

S = J dtL (52) 

without any metric factors in the integrand. Hence the 
Lagrangian for a point particle in a gravitomagnetic field 
is given by the square bracket of the matter term in 
Eq. jSU>- 

The canonical momentum of Lagrangian mechanics is 
defined by (p can )k — dL/dx k , where k = 1,2,3. From 
Eq. I|51|) we obtain in Cartesian coordinates 

Pcan = m(x + Ag). (53) 

This is the same equation as in classical Lagrangian me- 
chanics for point particles in an electromagnetic field, 
except that the electric charge q is replaced by the mass 
m. — The kinetic momentum is mx. It can be used as the 
input for solving Einstein's equations, because it is di- 
rectly determined by measurements without prior knowl- 
edge of the gravitomagnetic vector potential Ag = (3, 
which is an output of solving Einstein's equations. On 
the other hand the canonical momentum (for a given 
measured state of motion) depends on the gravitomag- 
netic vector potential A g . Therefore the canonical mo- 
mentum cannot be used as an input for solving Einstein's 
equations. The canonical momentum cannot be deter- 
mined by a FIDO from matter measurements without 
prior knowledge of the 4-geometry. — In curvilinear coor- 
dinates of Euclidean 3-space the general definition of the 
canonical momentum, (p ca . n )k = dL/dx k , gives 

ip can )k = mg^(x n +A^) k,n= 1,2,3. (54) 

From its definition via the Lagrangian, the canonical mo- 
mentum is a 1-form in 3-space, i.e. it has a lower 3- 
index. — On the other hand we can also start from the 
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4- velocity u", which is the archetype of a 4- vector (tan- 
gent vector), multiply it with the mass to obtain the 4- 

momentum p v , and pull down the 4- index with g£} . For 
the spatial components pk this gives pk = rn g^ju" — 

m 9^l(x n + j4") for v <C c, which is the same as the 
canonical momentum in Eq. I|54(l . — We note that in 
spherical coordinates p^ has the physical meaning of 
canonical angular momentum around the z-axis. 

Mathematically one can freely raise and lower indices 
of the vectors p resp. v using the metric tensor, or con- 
vert e.g. from lower indices to LONB components us- 
ing p^ = (ea)^P M . But physical quantities are prototypes 
for either LONB components, or 1-form components, or 
contravariant vector components. For exemplc in spher- 
ical coordinates the angular velocity around the z— axis 
is given by an upper-index component, = dx^/dt, 
the canonical angular momentum around the z— axis is 
given by a lower-index component, p^ = mg^^dx^ / dt = 
mr 2 sin 2 9 [(dx^ /dt) + A^], while the measured linear mo- 
mentum in the direction is given by a LONB com- 
ponent, p^ = m^jg^dx^ I dt = mr sin 9 (dx^ / dt) . The 
physics question dictates, which type of index is the rel- 
evant one. 

In the special case of azimuthal symmetry, i.e. when 

= £ is a rotational Killing vector, the canonical angu- 
lar momentum p^ — (p, £) is conserved, while the mea- 
sured, kinetic angular momentum of matter, (r sm9p ( j ) ), 
is not conserved because of the gravitoelectric induction 
field E g , Eq. I|29l) . The canonical angular momentum is 
relevant for the time-evolution, i.e. for the dynamics and 
conservation laws, not for kinematics (measurements at 
a given time). 

In the continuum description of matter (energy cur- 
rent density) the LONB components T^f. can be used as 
input for solving Einstein's equations, since they can be 
measured without knowing the output (3j. of Einstein's 
equations. On the other hand the coordinate-basis com- 
ponents T° k cannot be used as an input, because they 
cannot be determined by matter measurements without 
prior knowledge of Pk- 

In Mach's principle the input is the observed angular 
velocities of cosmic matter around us as in Eq. resp. 
the measured kinetic angular momenta (p+p)(rx v) as in 
Eq. I|42ll . This has also been recognized by Lynden-Bcll 
et al |l6(, who wrote (in section 2.4) "In Mach's principle, 
we wish to relate the dragging not merely to the angular 
momentum distribution but also to the distribution of 
angular velocity of the 'stars' about us." Unfortunately 
Lynden-Bell et al |l6j | did not follow up on this important 
observation, i.e. they did not solve the corresponding 
equation. Instead, in the rest of their paper, they used 
the canonical angular momentum T° as the input on the 
right-hand side of Einstein's equation. For the reasons 
explained in this section, measured kinetic matter input 
versus dynamical-geometric output of solving Einstein's 
equations, we think that their approach would be a very 
large detour to obtain our Eq. (J2J, which is the most 



direct expression of the dragging of inertial frames. — If 
one follows their proposal, the (H A E ) — term on the left- 
hand side of the Einstein equation Q37[) is absent, and 
there is no exponential suppression factor multiplying the 
canonical angular momentum distribution. This version 
of the momentum constraint equation may be may be 
relevant in other contexts. — In response to our confer- 
ence paper |l8| Bicak, Lynden-Bell, and Katz |20| gave 
arguments why they thought that the canonical angu- 
lar momentum density should be considered more fun- 
damental, although they agreed that "the problem with 
given angular velocities may be closer to Mach's original 
principle." They solved the equations for both input pos- 
sibilities (canonical versus kinetic angular momentum) 
for k = 0,±1. 



IX. EINSTEIN'S OBJECTION: "IF YOU HAVE 
A TENSOR T M „ AND NOT A METRIC ..." 

Einstein's objection to Mach's Principle in his letter to 
Felix Pirani of 2 February 1954 [32], as quoted by Ehlers 
in H3: 

"If you have a tensor and not a met- 
ric, then this does not meaningfully describe 
matter. There is no theory of physics so far, 
which can describe matter without already 
the metric as an ingredient of the description 
of matter. Therefore within existing theories 
the statement that the matter by itself deter- 
mines the metric is neither wrong nor false, 
but it is meaningless." 

From this argument Einstein drew the conclusion that 
"one should no lon ger speak of Mach's principle at all", 
quoted by Rcnn in 33] . 

We agree with Einstein's statement as quoted by 
Ehlers in |33| . as long as one refers to coordinate-basis 
components T^„. Einstein's valid criticism applies to T ^ 
in the proposal of Bicak et al [2(j ■ On the other hand, for 
LONB components T^ b our conclusions are entirely differ- 
ent. We have already explained this issue in the context 
of linear gravitomagnetism in the previous section. We 
shall now discuss Einstein's objection in a general con- 
text. 

The metric in LONB components, the Lorentz metric 
?7 a g = diag(— 1, +1, +1, +1), by itself does not encode spe- 
cific information about the particular Riemann space at 
hand (e.g. perturbed FRW metric versus Schwarzschild 
metric) , while the metric in coordinate-basis components 
gaflix) encodes all specific information about the partic- 
ular Ricmannian space at hand (curvature etc), i.e. it 
is metric data specific to the particular geometry. Both 
objects, the universal Lorentz metric rj &0 and the case- 
specific data-set g a p(x), are metrics, but they play totally 
different roles in the context of input versus output when 
solving Einstein's equations. 
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In Cartan's LONB formalism the case-specific data 
set which encodes the specific metric information (about 
the particular Riemann space at hand) is obtained by 
measuring the LONB components of the coordinate ba- 
sis vectors all over the particular Riemann space, 
i.e. measuring (e M )£, which gives g^vix) directly via 

9nv{x) = (e M )£ (e u ) h x r^g. 

In the LONB formalism a choice must be made for a 
specific field of LONBs, ea(P). For a coordinate system 
P — ► Xp given on a manifold one class of choices for 
LONBs is obtained by an orthonormalization procedure 
starting from the coordinate basis, as we have done in 
Eq. JT5J). 

Having the universal metric rj^ available is analogous 
to a surveyor having meter sticks available as a universal 
tool before starting the work of surveying a specific co- 
ordinatized two-dimensional Riemann surface. The case- 
specific output data-set of the surveyor's work is a table 
of measurements of (e^) p , hence g^ v {P) for a grid of 
points. 

The basic issue is not the fact that the universal metric 
tool r/«g makes it possible to perform the work of mea- 
suring the case-specific data-set g^(P) all over space- 
time. The basic issue is that such measurements are not 
needed as a case-specific input data-set, if one wants to 
solve Einstein's equations and thereby predict g fiu (x) as 
a case-specific output data-set. 

The crucial question is: What are the needed case- 
specific "ingredients" , what is the needed case-specific in- 
put data-set for solving Einstein's equations and thereby 
obtaining the case-specific metric output data-set? 

The input data set for solving Einstein's equations is 
obtained by performing the work of measuring the LONB 
components of the 4-momenta p of all particles in tan- 
gent spaces all over space-time, i.e. measuring p a . Sum- 
ming over all particles within one measured cm 3 around a 
point gives the LONB components T-g(P), the observed 
matter input data for solving Einstein's equations. To 
obtain this matter input one only needs the universal 
Lorentz metric 77^ = diag{— 1, +1, +1, +1}, but one does 
not need the case-specific geometric data g a p(x), which 
are the case-specific geometric output data from solving 
Einstein's equations and not available at the input level. 

We conclude that Einstein's objection quoted above 
does not apply to LONB components T.g, because no 
case-specific metric data-set is needed as an ingredient to 
obtain T a g. — If one wants to implement the phrase "mea- 
sured matter properties (by themselves) act on space, 
telling it how to curve " |24j. one is forced to use LONB 
components. And for LONBs the most convenient for- 
malim is the one of E. Cartan. 



X. LOCAL VORTICITY MEASURED BY 
NON-ROTATING OBSERVERS 

Local vorticity is defined as curlvfluid measured in the 
local inertial coordinate system which is comoving with 



the fluid, i.e. measured relative to the axes of local gy- 
rosopes ("local compass of inertia"). In general coordi- 
nates the local vorticity measured by non-rotating ob- 
servers is 

uj a = -e^upV-yUs, (55) 

where g 0123 = — 1. 

Mach's principle, formulated as a general hypothesis 
by Mach [5| [3(J] and made precise in Eq. states 
that the gyroscope axes here follow the rotational energy 
currents J e of matter in the universe averaged with a r~ 2 
weight and an exponential cutoff at the P-dot radius. In 
general the gyroscope axes here most definitely do not 
follow the motion of the local fluid here, i.e. relative 
to gyroscopes the local vorticity is nonzero according to 
Mach's principle in general. 

There is a special model, a rigid rotation of a fluid out 
to a rotation radius P ro t- In the limit P ro t/ 1 RdH/dt ~ * 00 
this special model produces an unperturbed FRW uni- 
verse, where the local vorticity measured by non-rotating 
observers vanishes. 

Unfortunately many authors have considered the van- 
ishing of the vorticity relative to the local compass of in- 
ertia to be a test for Mach's principle. See e.g. Ozsvath 
and Schucking's 01 solution and discussion of a Bianchi 
IX model for perturbations of the closed and static Ein- 
stein universe. In contrast we conclude that the vanishing 
of the vorticity relative to the local compass of inertia is 
not relevant as a test of Mach's principle. 

The claim of Ozsvath and Schiicking [3J| that their 
Bianchi IX model violates Mach's principle is incorrect, 
because their claim is based on the non-vanishing of the 
local vorticity. Now the question is: Does the Bianchi IX 
model of Ozsvath and Schiicking satisfy Mach's princi- 
ple (which involves an integral over energy flows in the 
universe)? 

XI. BOUNDARY CONDITIONS FOR 
EINSTEIN'S EQUATIONS IN ASYMPTOTIC 
MINKOWSKI SPACE 

For problems in asymptotic Minkowski space Ein- 
stein's equations by themselves are not sufficient to deter- 
mine a solution, they must be supplemented by boundary 
conditions at spatial infinity. 

An example is the Schwarzschild solution. Consider 
Earth in asymptotic Minkowski space. In General Rela- 
tivity we are free to choose any coordinate system. We 
choose spherical coordinates with Earth-fixed axes, but 
all the same we write down the standard form of the 
Schwarzschild solution. This certainly satisfies Einstein's 
differential equations for the Earth as a source. But it 
is not the correct solution. Just consider a geostation- 
ary satellite. With this solution the equation of motion 
(geodesic equation) predicts that the satellite will fall 
down. What is missing are the centrifugal forces, which 
arise from our choice of Earth- fixed axes. 
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We now can do one of two things: 1) Either we do not 
admit coordinate systems which (in the asymptotic re- 
gion for the solar system) are rotating relative to quasars. 
Then the usual asymptotic condition g^ v — > rj^ v is adae- 
quate. But excluding a class of coordinate systems is 
against the spirit of general relativity. 2) Or we do admit 
coordinate systems which are rotating asymptotically rel- 
ative to quasars, e.g. we admit coordinate systems with 
Earth-fixed axes. But then we need encoding boundary 
conditions in the asymptotic part of the solar system. 
The boundary conditions encode fictitious forces (e.g. the 
Coriolis force, which is equivalent to a homogeneous grav- 
itomagnetic field) in the asymptotic part of the solar sys- 
tem. In other words the boundary conditions encode the 
effects of sources outside the system considered (outside 
the solar system), they encode the influence of the cos- 
mological sources. 

The necessity of boundary conditions at spatial infinity 
is familiar from electromagnetism. Boundary conditions 
arc needed for fields, whose sources lie outside the system 
under discussion. In ordinary magnetism the sources ex- 
ternal to the system could be coils with electric currents 
generating applied external magnetic fields for the exper- 
iment under consideration. 

What is the general form of the boundary conditions in 
the region asymptotic to the solar system? The forms of 
the asymptotic force fields are familiar from the fictitious 
forces in classical mechanics arising in coordinate systems 
which are rotating relative to inertial systems. In the 
vector sector of linear perturbations (gravitomagnetism) 
two fictious forces arise, since the centrifugal force m[0 x 
[fx f2]] is absent at the level of linear perturbation theory: 

1) The Coriolis force 2m[v x f2], which corresponds to 
a homogeneous gravitomagnetic field, 

B g (2.*)">(4)EL = -2fi$o, (56) 

where 1S the precession vector of the spin axes of 

asymptotic gyroscopes (i.e. axes of inertial frames) rela- 
tive to the axes of our FIDO's (e.g. Earth-fixed axes). 

2) The fictitious force due to a non-constant angular 
velocity, m[f x dQ,/dt], which corresponds to a gravito- 
electric induction field, 

S g (2,t)^{S g )"" = [?x±fi%y o ]. (57) 

The asymptotic B g and E g fields are solutions of the 
source-free Einstein equations, i.e. they do not have any 
sources within our system (solar system). Mach's princi- 
ple implies that the sources for these fields are given by 
cosmological matter (which is at infinity from the point of 
view of the solar system), as we have shown in Eqs. (|42l 

El- 
One must distinguish two types of boundary condi- 
tions: 1) Boundary conditions which allow to include 
and encode the contribution of sources outside the sys- 
tem considered. 2) Boundary conditions which exclude 



unphysical solutions, e.g. the exponentially growing so- 
lution of Eq. lptU|) . Boundary conditions discussed in the 
context of Mach's principle refer to the encoding type. 

In this paper we have shown that in cosmological grav- 
itomagnetism there is no need for boundary conditions 
of the encoding type, because of the exponential cutoff 
for the contributions by the measured energy-momentum 
tensor of matter. 

Einstein and others had hoped for some time that one 
could implement the ideas of Mach by specifying some 
boundary conditions on th e eq uations of General Rela- 
tivity. But in 1934 Einstein [35| concluded that the prob- 
lem of boundary conditions (for Mach's principle) could 
only be solved by going to a spatially closed universe: 

"In my opinion the general theory of rela- 
tivity can only solve this problem [of inertia] 
satisfactorily if it regards the world as spa- 
tially self-enclosed." 

In contrast we have shown in this paper that General Rel- 
ativity for linear cosmological perturbations of a spatially 
flat FRW universe satisfies Mach's principle without the 
need for boundary conditions of the encoding type. In a 
companion paper |23l ] we shall show that this is also true 
for perturbations of FRW universes with k = — 1, i.e. for 
open, spatially hyperbolic universes. 

XII. MACH'S PRINCIPLE IN THE PRESENCE 
OF SCALAR AND TENSOR PERTURBATIONS 

Linear scalar perturbations cause weak gravitational 
lensing. The geodesies from the gyroscope to quasars (on 
St) will suffer geodesic deviation. This forces us to adapt 
the method of fixing the spatial axes of our FIDO's, i.e. 
our LONB's, as follows: We demand that the spatial axes 
of the LONB's do not rotate relative to the average of the 
directions of geodesies (on E t ) to distant quasars, where 
the weight function in the average is the toroidal vector 
field X7_ 1 = —iUYii—x^n |23l |. which is tangential to 
the 2-sphere, 

= J dr W(r) J dn (X~ =l m y ■ v^. (58) 

The FIDO at the origin determines v^y 5 from the mea- 
sured angular velocities of tangents to geodesies on Sj 
from the origin to quasars in the asymptotic FRW uni- 
verse. An observational window function in the asymp- 
totic FRW universe is denoted by W(r). 

The vector field of distortions by weak gravitational 
lensing belongs to the scalar sector. This vector field is 
orthogonal (under integration over the 2-sphere) to the 
toroidal vector spherical harmonics X^ =l m , because the 
two fields belong to the opposite parity sequence. The 
scalar sector has P — (— 1) , the toroidal vector spherical 
harmonics have P = (— . Therefore scalar pertur- 
bations cannot affect Mach's principle, if the gyroscope 
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axes are measured relative to the average of the direc- 
tions of geodesies to asymptotic quasars with the weight 
function X^ l m . 

Gravitational waves at the linear level in the perturbed 
region between asymptotic FRW space and the gyroscope 
also cannot affect the average of geodesies in Eq. (l5%ll . 
again because the tensor sector is orthogonal to the vec- 
tor sector. 

Gravitational waves going through the gyroscope can- 
not produce a gyroscope precession (relative to quasars 
in the asymptotic unperturbed universe). This is be- 
cause the precession of the gyroscope only depends on 
the three independent rotation coefficients (w«)g for a 
displacement of unit measured time along the worldline 
of the gyroscope. These three rotation coefficients form 
an axial 3-vector, which cannot be affected by a 3-tensor 
(gravitational wave sector) for reasons of symmetry un- 
der 3- rotations. 

If angular velocities are measured relative to the gy- 
roscope axes (not relative to distant quasars), one has 
the value figyro — on the left-hand side of Eq. (|42[1 . 
and there is no need to modify the formula due to scalar 
or tensor perturbations, because Eq. (|42(l automatically 
performs a projection on {£ = 1, odd parity sequence}. 

Einstein in his paper "Prinzipicllcs zur allgcmeinen 
Relativitatstheorie" of 1918 coined the term "Mach's 
Principle" and gave the following formulation for it: 

"The G- field is entirely determined ... by 
the energy tensor of matter." 

The G- field refers to the metric tensor field g fiv (x). 

For linear cosmological perturbations Einstein's for- 
mulation is too strong to be valid, and unnecessarily 
strong for Mach's original purpose of explaining (giving 
the physical cause for) the evolution of groscope axes 
(i.e. non-rotating frames): 1) Too strong, as has been 
stated many times, because g^ v can contain gravitational 
waves, which can exist without matter being present. 
2) Unnecessarily strong, as we shall now explain: Gravi- 
tational waves cannot affect gyroscope axes directly, be- 
cause a tensor perturbation cannot affect the spin evolu- 
tion, dS/dt, which is an axial vector. And gravitational 
waves cannot affect the {X^ =1 m )-average of the geodesies 
from a gyroscope to the asymptotic unperturbed quasars. 

In our approximation of linear cosmological pertur- 
bations a third question remains unaddressed: Since 



the gyroscope's precession is determined by the B g 
field at the gyroscope's position, is B g determined by 
an equation which has a gravitational pseudo-angular- 
momentum tensor (in addition to the angular momentum 
of matter) as a source term? Such an equation could not 
be Einstein's equation in the standard way of separating 
output curvature terms from the matter input. 

XIII. OUTLOOK 

1. The analysis of this paper has been extended to 
FRW universes with k = +1 (spatially closed, 
spherical) and with k = — 1 (spatially open, hy- 
perbolic) with analogous results. This will be pre- 
sented in a companion paper p3| . 

2. Causality and astronomical observations make a 
formulation desirable, where B g at the space-time 
point P is given by an integral over sources on the 
backward light-cone of P. 

3. Second order perturbation theory: Does perfect 
dragging of gyroscope axes still hold? What is 
the contribution of gravitational waves to dragging 
of gyroscope axes? — Can one derive exact results 
about perfect dragging of gyroscope axes? 

4. Does the Bianchi IX model of Ozsvath and 
Schiicking satisfy Mach's principle (which involves 
an integral over energy flows in the universe)? 

5. Accelerated frames of reference include not only ro- 
tating frames (extensively discussed by Mach), but 
also frames with linear acceleration relative to iner- 
tial frames. Docs cosmological matter (with an ex- 
ponential cutoff) cause perfect frame dragging with 
respect to linear accelerations? 
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